clear all

set more off
set scheme modern

if c(username)=="" { //insert username
cd "" // insert root path
}

if c(username)=="jmc4qg" {
cd "/Users/jmc4qg/The Lab Dropbox/Jonathan Colmer/Intergenerational_Presentation/Submissions/JPE - Micro/Figure_Replication/"
}

** Fig 1a **

import delimited "Data/Fig1.csv", encoding(ISO-8859-1) clear
keep if v7 == "Fig1a.csv"
drop v1 v4 v5 v7

rename v6 year
rename v2 estimate
rename v3 stderror

sort year

gen u = estimate + 1.96*stderror
gen d = estimate - 1.96*stderror

twoway(line u year, fcolor(black) lcolor(black) lw(thin) lpattern(dash))(line d year, fcolor(black) lcolor(black) lw(thin) lpattern(dash))(line estimate year, lcolor(black) lpattern(solid) lwidth(medthick) legend(off) ///
ytitle("Change in TSP ({&mu}g/m{sup:3})", size(large)) xtitle("Year", size(large)) graphregion(color(white)) plotregion(color(white)) xlabel(1969(1)1975, labsize(large)) xline(1971) ///
ylabel(-30(10)10, labsize(large)) yline(0))
graph export "Figures/Fig1a.pdf", replace


** Fig 1b **

import delimited "Data/Fig1.csv", encoding(ISO-8859-1) clear
keep if v7 == "Fig1b.csv"
drop v1 v4 v5 v7

rename v6 year
rename v2 estimate
rename v3 stderror

sort year

gen u = estimate + 1.96*stderror
gen d = estimate - 1.96*stderror

twoway(line u year, fcolor(black) lcolor(black) lw(thin) lpattern(dash))(line d year, fcolor(black) lcolor(black) lw(thin) lpattern(dash))(line estimate year, lcolor(black) lpattern(solid) lwidth(medthick) legend(off) ///
ytitle("Change in TSP ({&mu}g/m{sup:3})", size(large)) xtitle("Year", size(large)) graphregion(color(white)) plotregion(color(white)) xlabel(1969(1)1975, labsize(large)) xline(1971) ///
ylabel(-30(10)10, labsize(large)) yline(0))
graph export "Figures/Fig1b.pdf", replace

** Fig 1c **

import delimited "Data/Fig1.csv", encoding(ISO-8859-1) clear
keep if v7 == "Fig1c.csv"
drop v1 v4 v5 v7

rename v6 year
rename v2 estimate
rename v3 stderror

sort year

gen u = estimate + 1.96*stderror
gen d = estimate - 1.96*stderror

twoway(line u year, fcolor(black) lcolor(black) lw(thin) lpattern(dash))(line d year, fcolor(black) lcolor(black) lw(thin) lpattern(dash))(line estimate year, lcolor(black) lpattern(solid) lwidth(medthick) legend(off) ///
ytitle("Change in TSP ({&mu}g/m{sup:3})", size(large)) xtitle("Year", size(large)) graphregion(color(white)) plotregion(color(white)) xlabel(1969(1)1975, labsize(large)) xline(1971) ///
ylabel(-30(10)10, labsize(large)) yline(0))
graph export "Figures/Fig1c.pdf", replace

** Fig 1d **

import delimited "Data/Fig1.csv", encoding(ISO-8859-1) clear
keep if v7 == "Fig1d.csv"
drop v1 v4 v5 v7

rename v6 year
rename v2 estimate
rename v3 stderror

sort year

gen u = estimate + 1.96*stderror
gen d = estimate - 1.96*stderror

twoway(line u year, fcolor(black) lcolor(black) lw(thin) lpattern(dash))(line d year, fcolor(black) lcolor(black) lw(thin) lpattern(dash))(line estimate year, lcolor(black) lpattern(solid) lwidth(medthick) legend(off) ///
ytitle("Change in TSP ({&mu}g/m{sup:3})", size(large)) xtitle("Year", size(large)) graphregion(color(white)) plotregion(color(white)) xlabel(1969(1)1975, labsize(large)) xline(1971) ///
ylabel(-30(10)10, labsize(large)) yline(0))
graph export "Figures/Fig1d.pdf", replace
